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We present a derivation of the low energy Lagrangian governing the dynamics of the spin degrees 

of freedom in a spinor Boso condensate, for any phase in which the average magnetization vanishes. 
This includes all phases found within moan-field treatments except for the ferromagnet, for which 
the low energy dynamics has been discussed previously. The Lagrangian takes the form of a sigma 
model for the rotation matrix describing the local orientation of the spin state of the gas. 



I. INTRODUCTION 
A. Spin ordering in ultracold gases 

In recent years the field of ultracold atomic physics 
has attracted the attention of a great many condensed 
matter theorists, largely due to the prospect of finding 
novel realizations of many-body systems. Part of the ap- 
peal doubtless lies in the exquisite experimental control 
that may be exercised over the parameters of a system in 
which many of the complicating factors familiar from the 
solid state (disorder, phonons, etc.) are absent. Certain 
intrinsic aspects of ultracold atomic gases - not depen- 
dent on the spec;ifics of the experimental setting - are 
qualitatively novel, however, and without antecedent in 
the study of condensed matter. In this latter category we 
may place the possibility of spontaneous ordering of the 
spin degrees of freedom in a Bose gas. Indeed, prior to the 
'ultracold revolution' the only Bose superfluid that could 
be studied in the laboratory was ^He, which has zero spin. 
With the advent of optical trapping of Bose condensates 
of alkali atoms, which allows for a fully rotationally in- 
variant setting, the experimental study of spin ordering 
within a hyperfine multiplet came within reach^'^. 

The earliest theoretical works motivated by these de- 
velopments explored possible ordered phases using a 
mean field description, in both the spin 1 and spin 2 hy- 
perfine multiplets^"^ . In this description the concept of 
spontaneous symmetry breaking plays a central role. Up 
until very recently, however, there were no experiments 
in which this spontaneous ordering was apparent. The 
reason for this is that the simplest experimental protocol 
for the investigation of the hyperfine state of the gas is to 
apply a magnetic field gradient to split a gas cloud into 
different components in a Stern-Gerlach experiment 
The different components are subsequently imaged to de- 
termine their (relative) occupancies. This technique nat- 
urally imposes a quantization axis, and any information 
concerning the coherence between different hyperfine lev- 
els is lost. Thus magnetic alignment in the plane perpen- 
dicular to this axis, which depends on the relative phase 
of these different levels, cannot be observed. 

The characterization of magnetic ordering in atomic 
gases has taken a leap forward in the last few years thanks 
to the work of the Berkeley group, who demonstrated in 
situ dispersive imaging of the transverse magnetization of 



a gas of ^''Rb in the spin 1 multiplet^, and subsequently 
employed this technique to investigate a number of fas- 
cinating aspects of this system, including the dynamics 
of spontanteous symmetry breaking, defect production, 
and the role played by magnetic dipole forces^" "'^'^ . 

The above developments illustrate two important 
needs. Firstly, imaging of the spin order was necessary 
to bring much of this new physics out into the open. Sec- 
ondly, the nonequilibrium character of most experiments 
reriuires that the mean field theory of equilibrium ordered 
states be supplemented with a dynamical description of 
the relevant order parameters. 

It is our hope that the next few years will see the devel- 
opment of imaging techniques capable of detecting some 
of the spin orders to be discussed in Section III, of which 
an average magnetization is only the simplest. The aim 
of this work is to address the second need: uncovering 
the low energy description of the order parameter. For 
the case of the Bose ferromagnet, which is appropriate to 
the spin 1 ^^Rb system, this description was provided in 
an earlier paper^^. In this work we will focus instead on 
states with vanishing average magnetization. For reasons 
that will be become clear in the following sections these 
cases are qualitatively different. 



B. Low energy descriptions 

In an ordered phase of matter we expect that the low 

energy degrees of freedom consist of variations of the or- 
der parameter on some manifold of symmetry broken 
states (Goldstone modes), together with any conserved 
quantities. In our earlier work on the Bose ferromagnet^^ 
the degrees of freedom were the local magnetic moment 
m and the superfluid velocity v. In the long wavelength 
limit these were found to obey the coupled equations 
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-^ + v'V is the convective derivative. Eq. (1) 
is a modified Landau-Lifshitz equation and gives rise 
to quadratically dispersing spinwaves when linearized on 
around a solution m = const.. 
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The quadratic dispersion is a consequence of the two 
transverse deviations of the order parameter being canon- 
ically conjugate. For the phases with vanishing average 
magnetization that are the focus of this work, the conju- 
gate variables involve deviations from the order param- 
eter manifold. This results in linearly dispersing Gold- 
stone modes. The situation is analogous to the case of 
spin waves in an antiferromagnet, where the conjugate 
variables are the difference in magnetization on neigh- 
boring sites - the Neel order parameter - and the sum. 

It also follows from the vanishing of the average magne- 
tization that the order parameter dynamics and super- 
fluid flow are decoupled except for a global topological 
constraint. Again, this is quite different from the case of 
the ferromagnet, where the two are coupled together in 
the equations of motion. It follows that we can write a 
Lagrangian for the spin degrees of freedom only. This La- 
grangian is expressed in terms of a rotation matrix R that 
specifies the local orientation of the spin state relative to 
some reference state. By expressing the matrix elements 
of R in terms of an orthonormal triad Rab = {^b)a' with 
Ca ■ = Sab, the spin Lagrangian may be written 
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Here are Qa are some constants to be specified later, 
which depend on the ordered phase in question. The La- 
grangian Eq. (3) is the main result of this paper. The 
only other phase for which the spin Lagrangian was pre- 
viously obtained is the polar phase of the spin 1 gas, to 
be discussed below, in which = = 0, a = 2, 3, giving 
the familiar 0(3) model^^. 



C. Outline of this paper 

The structure of the remainder of this paper is as fol- 
lows. In Section II we review the basic description of 
a spinor condensate, and in particular the structure of 
the interaction Hamiltonian. Section III describes in de- 
tail the ground state manifolds of the different phases 
of spinor condensates, beginning with the simplest non- 
trivial case, spin 1, before introducing the Majorana (or 
stellar) representation that is very useful in visualizing 
spin states. The global structure and local geometry 
of the order parameter manifolds are then introduced, 
as well as a parametrization for the dynamical variables 
conjugate to the order parameter. After this the deriva- 
tion of the low energy Lagrangian is a relatively simple 
matter, and is described in Section IV along with the 
derivation of the equations of motion. We have tried to 
keep the presentation pedagogical throughout, though at 
various points there are parenthetical technical comments 
that the casual reader should feel free to ignore. 

In a related work, Barnett et al}^ have derived the full 
equations of motion of a spinor condensate in terms of the 
Majorana representation and applied a group-theoretical 



analysis to the determination of all normal modes about 
an ordered ground state. We pursue a complimentary 
goal of obtaining the full nonlinear Lagrangian for the 
Goldstone modes only. 



II. BASICS OF SPINOR ROSE CONDENSATES 
A. Lagrangian 

One could take the point of view that the descrip- 
tion of the dynamics of a dilute spinor gas is no dif- 
ferent from its spinless counterpart, being governed by 
the time-dependent Gross-Pitaevskii equation, with La- 
grangian density 



(4) 



(we set = m = 1 from now on). For a spin s gas 4> is 
a 2s -h 1 component spinor and the Hamiltonian density 
has the form 
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where the first term is the kinetic energy. The interaction 
part 7iint(0^, </*) is quartic in ({> and its form will be given 
below for s = 1 and 2. We will not discuss the influence 
of the trapping potential, save to assume that it preserves 
rotational symmetry. For stationary solutions of the form 
(f){r,t) = e^*^*(/)(r) the time-dependent description re- 
duces to the time-independent Gross-Pitaevskii theory, 
with /i the chemical potential. Instead of treating the 
action S = J dvdt C classically, we can interpret it as the 
quantum action in a coherent state path integral. Little 
that we will have to say will depend upon this distinction. 

In fact this superficial similarity between the spinless 
and spinful problems is quite misleading. The ground 
state in a uniform system corresponds to some constant 
(j). In the spinless case this fixes ^ up to a phase once 
the density p = (f)^ (f) is specified. But in the spinful case 
we still have to find the correct 'direction' of (p in the 
complex 2s -|- 1 dimensional spinor space, determined by 
minimizing the interaction Hamiltonian Wint(<?!'^5 <i>)- The 
interactions will be assumed to respect rotational sym- 
metry, so this minimum is only unique up to rotation 
(specified by the three Euler angles, say). Choosing this 
rotation - starting from some arbitrary reference state 
- specifies the spontaneous breaking of rotational sym- 
metry in the ground state. The slow variation of this 
rotation in space and time constitutes the low energy 
dynamics of the system, the description of which is the 
subject of this work. To characterize these low energy 
manifolds we must first specify the form of Hint- 



R. Interaction Hamiltonian 

The structure of Tiint has been discussed in many 
works, starting with the first papers treating spinor Bose 
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condensates. We will therefore keep the following dis- 
cussion relatively brief. Low energy scattering between 
a pair of bosons occurs in the s-wave channel only, and 
can be treated as a ^-function interaction, characterized 
by a set of interaction constants gs, 5 = 0, 2, ... 2s for a 
pair of bosons with total spin S. Bose symmetry dictates 
that the interaction vanishes for odd total spin. 

For spin 1, the resulting interaction may be presented 
in the form^'"* 

Hint = I (0^0)% 1(0^5^0)' (6) 

where S''^*^ are the spin s angular momentum matrices. 
Eq. (6) is the sum of a density-density and a spin-spin 
interaction. For spin 2 we have^'^^ 

Here cj) ■ (j) = Yl,m=-s (^1)"^™ '/'m'/'-m is a scalar repre- 
senting the amplitude of singlet pairs of spin s (Eq. (6) 
can be expressed using this operation instead of the 
spin-spin interaction. In the spin 2 case both terms 
are needed). Our use of C2 for two different quantities 
matches the notation of the works cited above, where one 
may find explicit expressions for Cj i = 0, 1, 2 in terms the 
constants gs defined above. 

In the spin 1 case it is fairly clear how to minimize 
Eq. (6) at fixed density. For the spin 2 case things are less 
obvious. In the next section we will discuss a method of 
parametrizing the spinor (j) to make this operation trans- 
parent. 



III. IDENTIFYING THE SPIN DEGREES OF 
FREEDOM 

In this section we will introduce the Majorana (or stel- 
lar) representation of spin states. This provides a vivid 
way to picture spin ordering in higher spin condensates in 
which rotational symmetry is manifest. None of the cal- 
culations of Section IV depend upon this representation; 
its use is rather in providing a concrete way to picture 
the ground state manifold. 

Before beginning it is worth setting the scene with a 
more pedestrian discussion of the spin 1 case^'^. 



A. Phases of the spin 1 gas 

Let us minimize Eq. (6) with a spinor (j) normalized to 
unity (thus we are adopting units in which the density 
p = 1). This is a matter of minimizing (maximizing) 
(j)^S'^^^(j> for C2 > (c2 < 0). One way to make the 
resulting states more clear is to write the complex spinor 
(j) = a + ib where a and b are two real vectors satisfying 
+ b"^ = 1, and to work in cartesian components. The 



relationship to the usual components (j)m m = —1, 0, 1 is 

(^^x = (^1 - (p-i) 

"Aj/ = {<t>i + 4>-i) 

<t>z = <i>o (8) 

In this basis the angular momentum matrices take the 
form (S'l^M = — iejjfe.Then we have 

(t>^S'^^U = 2a X b 

For C2 < the interaction energy is minimized for a 
and b perpendicular and equal in magnitude. This state 
is termed ferromagnetic as it corresponds to maximal 
polarization of the spin. The resulting order parameter 
manifold c;orresponds to the set of all configuration of a 
pair of orthogonal vectors, and is thus identified with the 
group of rotations 50(3) 

For C2 > a and b are aligned. The resulting polar 
state (the name originates from an analogous state in 
superfluid ^He) can therefore be written as 

<Ar^,e = e'^n (9) 

for n a real unit vector. Note that this paranietrization 
has some redundancy in that the points [9 + tt, —n) and 
{0, n) are identified. The resulting manifold is known as 
the mapping torus of the antipodal map of the sphere S'^. 

The global topology of the above ferromagnetic and 
polar order parameter manifolds naturally determines 
the character of the topological defects in the ordered 
phases, and certain features of the ordering transitions, 
some of which have already been discussed in the liter- 
ature''^^. This is not the focus of the present work and 
topology will not be further discTissed, even though the 
defect physics of the higher spin condensates promises to 
be highly non-trivial 

The polar state has e'^^^Vn e = 0- Nevertheless 
the above discussion makes it clear that polar ordering 
involves a choice of axis: the spinor (t>nfi is the m = 
state with respect to the axis n. It is natural to ask for 
an operator that acquires a non-zero expectation value 
in the polar state. The obvious candidate is the spin 2 
quadrupole tensor (or nematicity) 

= \ (S^'S^' + Si-'si^^) - '-^S^, (10) 

with 

Such expressions are familiar in the study of nematic liq- 
uid crystals, where the vector n is known as the director. 
In the liquid crystal context the identification of n and 
—n without the phase factor in Eq. (9) (the order pa- 
rameter manifold is then the real projective plane MP^) 
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FIG. 1: (color online) Majorana representation of a s = 3 
spin state, with arrows representing the principal spinors. 



makes for very different defect physics, however Ne- 
matic ordering in soHd state magnetic systems has been 
the subject of much experimental and theoretical work 
in recent years, with a good deal of uncertainty still re- 
maining. The observation of the polar state in the spin 
1 Bose gas would therefore be an important milestone. 

Searching for higher spin order parameters as the spin 
of the gas particles increases becomes arduous. We now 
turn to a more convenient representation of the spin or- 
der. 



B. Majorana (stellar) representation 



the corresponding 2s -I- 1 component spinor <j) is 

1/2 
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Next contract every index of ^ab - c with C^^ 



-s, . . . s. 



Denoting by ^(z) the resulting polynomial of order 2s, 
the fundamental theorem of algebra tells us 
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with M some normalization. Thus '^ab - c niay be writ- 
ten 



AB- 



■c = "(A^s • • • 7c) 



(12) 



with the principal spinors a a, /3b i etc. related to the 
{zi\ by 
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relations that are unchanged if we normalize the spinors, 
in which case they correspond to 2s points on the Bloch 
sphere with coordinates {0^, (^i}, i — 1, ... 2s 
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e*^^/2cos^,e-*^^/2sin^ 



The representation of a general spin s state that 
(sometimes) bears his name was discovered by Majorana 
in 1932^^, and independently several times since^'^'^"'^, 
though it has antecedents in 19**^ century mathematics^^. 
A very nice discussion can be found in Ref. 23. 

The result is very simple to state, and represents a gen- 
eralization of the Bloch sphere for spin 1/2 to arbitrary 
spin. Up to normalization and a phase - thus in more 
mathematical terms we are parametrizing the complex 
projective space CP^'' - an arbitrary spin s state can be 
specified by locating 2s indistinguishable points on the 
unit sphere (see Fig. 1). Such a configuration is some- 
times called a constellation, for reasons that will become 
clear. 

There are two steps to understand why this is so. 
First, imagine forming our spin s from 2s spin 1/2 in 
the totally symmetric subspace. An arbitrary state may 
then be written as a totally symmetric spinor <^ab - c — 
'^(AB - C)) where the round brackets denote the operation 
of symmetrization and each of the 2s -I- 1 indices can take 
the value t or ],. The relationship between <^ab - c and 



(7T.7i) = 



e^^-/2cos^,e- 



^^-/^sin^ 



(notice that <^ab - c in Eq. (12) is not in general normal- 
ized when the principal spinors are). Then (minus) the 
roots of ^{z) can be written Zi = e*'''' cot 9i/2 and corre- 
spond to stereographic projection from the north pole to 
the plane tangent to the sphere at the south pole. 

The beautiful feature of the Majorana representation 
is that rotations act simply as rotations on the Bloch 
sphere. Of course, it is useful to have an explicit ex- 
pression for the polynomial '^{z) in terms of the 2s -|- 1 
components of the spin s state 4>. It is easy to show 
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If spinor indices are raised and lowered using the an- 
tisymmetric tensors ^ab and (with e|| = 1 and 



^AB _ , \ 
e — —(-AB) 
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then rotational invariance upon contraction of indices is 
guaranteed (in fact the resuh is invariant under the larger 
group SL{2, C), a result that will be useful later) . If we 
denote by the result on cp of raising all indices of the 
corresponding symmetric spinor, then one can readily see 
that 4>m = (— l)*'''"*0-m5 and thus 

m=—8 

After raising indices of the principal spinors we have 
/a^ = —l/zi, etc.. Under stereographic projection 
z —1/z* represents the antipodal map on the unit 
sphere. Thus we see that the spinor ^* is represented by 
a set of points antipodal to those representing (p. Fur- 
thermore, the Majorana representation of a normalized 
state with |0 • (/>| = 1, corresponding to ^ = e^^^ip*, con- 
sists of pairs of antipodal points (and is thus only possible 
for integer spin). This fact will be useful in minimizing 
the interaction energy (recall the form of Eq. (7)). The 
transformation T : (pm {T4')m = 4>m in fact the 
(anti-unitary) operation of time reversal. 



1. Spin 1 

As mentioned in Section II B, Eq. (6), the interaction 
Hamiltonian in the spin 1 case may be written 

Win.=^(^vr-li0-0r (15) 

I 1 2 

For C2 > we should maximize . Based on the 

discussion of the previous section, this corresponds to 
placing the two points antipodally in the Majorana rep- 
resentation. It is evident that the corresponding spin 1 
spinor is just the symmetric m = state with respect to 
the resulting axis. This is just the polar state described 
before. 

I 1 2 

For C2 < I (/)•(/) I can be set to zero by making the two 
principal spinors equal since aACt^ = 0. It is clear that 

this represents the ferromagnet, being a maximally po- 
larized (or coherent) spin state, a result that generalizes 
to arbitrary s. 



C. Spin ordering in the Majorana representation 

The use of the Majorana representation to visualize 
spin ordering in a Bose gas was suggested in Ref. 24. Let 
us first see how the phases of the spin 1 gas discussed in 
Section III A appear in this representation, before moving 
on to the spin 2 case. 



2. Spin 2 

Suppose that in Eq. (7) ci > 0. The corresponding 
term in the Hamiltonian can be fully satisfied by states 
with = 0. By inspecting the spin 2 angular 

momentum matrices 
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it is easy to convince oneself that up to rotation the most 
general state satisfying these conditions is 
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(17) 



The principal spinors form a polyhedron with four iden- 
tical triangular faces known as a disphenoid (see Fig. 2). 
To fix the parameters 77, x in Eq. (17) we turn to the third 
term in Eq. (7). If C2 < 0, this term is minimized by plac- 
ing four points in two antipodal pairs so that |0- (/ip = 1. 



The principal spinors form a rectangle, corresponding to 
X = in Eq. (17). The aspect ratio of the rectangle 
varies with rj, with 77 = (2n -|- 1) 7r/3 being a square and 
rj = corresponding to a pair of points at either pole. 
The fact that the energy is minimized for any rj in this 
parameter regime is rather surprising, and we will return 
to it briefly below. This state will be referred to as rect- 
angular in the following (in Ref. 24 it was called nematic 
since jVa(, ^ 0). 

For C2 > |(^ • (^p should be minimized. This can be 
done by taking x = = 7r/2. The result is a regular 
tetrahedron. Note that this state has Mat = 0. 

For ci < one must compare the energy of the rectan- 
gular state with that of the ferromagnet. The resulting 
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FIG. 2: (color online) Top: disphenoid corresponding to the 
Majorana representation of the state with x — 7r/4, rj — n/2 
in Eq. (17). Bottom left: rectangular state with x = 0, r; = 
tt/2. Bottom right: tetrahedron with X ~ V ~ 7r/2. 

phase diagram can be found in Refs. 5,24. Note that 
while the rectangular state maximizes the magnitude of 
the quadratic scalar 

the tetrahedral state maximizes the cubic invariant 

J^ci>;^gci>g^ci>fl. (19) 

The inclusion of a sextic term proportional to | Jp in the 
interaction energy would lift the accidental degeneracy in 
the r] parameter discussed above. A microscopic deriva- 
tion of such a term is discussed in Refs. 25,26, with a 
positive sign favoring the square state [r] = {2n + 1) 7r/3) 
and a negative sign the uniaxial state in which two pairs 
of points coincide (?7 = 0) . The inclusion of such a term 
in our formalism is a straightforward matter, and we will 
not discuss it further. 

The rectangular and tetrahedral states evidently have 
certain discrete symmetries that are rather hard to dis- 
cern by inspection of Eq. (17), and indeed went unno- 
ticed in the earliest works on the spin 2 condensate^'^. 
It appears that the term 'cycHc' used in several works 
to describe the tetrahedral phase is a consequence of a 
misidentification of the symmetry. This illustrates the 
utility of the Majorana representation in the visualiza- 
tion of spin order. 



D. Ground state manifolds 

1. Global structure 

In Section III A we identified the order parameter man- 
ifolds of the phases of the spin 1 gas. With the help of 
the Majorana representation we can now do the same for 
the spin 2 case. Roughly speaking, we expect the man- 
ifold to consist of all configurations related to those of 
Section III C 2 by rotation. Some rotations will leave the 
configuration of principal spinors unchanged, however, so 
the manifold cannot simply be identified with the rota- 
tion group SO{3). In mathematical terms the problem 
is to determine the orbits of a reference spinor under the 
action of the spin s representation of the rotation group. 

If we ignore the phase of the spinor for a moment, so 
that we are considering orbits in CP^**, this problem can 
be solved using the Majorana representation by consid- 
ering the manifold of constellations generated by all pos- 
sible rotations^*'. For instance, the ferromagnetic spin 
1 state has orbit in CP^ equal to (and this is true 
in general for any spin s coherent state), while for the 
polar state we have with antipodal points identified: 
the real projective plane MP^. In general if one finds a 
configuration of principal spinors unchanged under some 
subgroup of r C 5*0(3) (the stabilizer subgroup of a con- 
stellation), then the orbit is given by S'0(3)/r. Thus in 
the case of spin 1, F = SO{2) for the ferromagnet and 
F — 0(2) in the polar phase, since in the latter case a 
parity transformation also leaves the points unchanged. 

The reader will notice that these are not the order 
parameter manifolds identified in Section III A for the 
spin 1 case. We have neglected the phase of the spinor, 
which is a real degree of freedom. We might then guess 
that any spinor on the order parameter manifold can be 
written 

^R,e = e''D^''>{R)cbo (20) 

where Z?(^^(i?) is the spin s representation of the rotation 
R and (j)Q is some reference spinor corresponding to the 
phase in question. This does not mean that the order pa- 
rameter manifold is 50(3) x U{1), because the stabilizer 
subgroups mentioned above leave the spinor correspond- 
ing to a particular constellation unchanged up to a phase. 
We denote these phases as A(7) , with 7 S F . They must 
form a one-dimensional unitary representation of F 

A(7i72) = A(7i) + A(72) mod 2tt. 

and allow us to make the identification 

4>R,e = 0fl7,0-A(7) 7 S F, (21) 

showing that the order parameter manifold is ^ 
where the tilde is to denote the action of F on 5*0(3) x 
U{1): (i?,0)^(i?7,0-A(7)) 

The simplest case to consider is the polar phase with 
F = 0(2), for which the only non-trivial phase is a —1 as- 
sociated with the parity transformation. If F were 50(2) 
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FIG. 3: (color online) Tetrahedron with 3-fold axis marked, 
see Eq. (23). 



with only trivial phases, the ground state manifold would 
be X S*^, with the first factor coming from the C/(l) 
and the second from S0{3) / SO{2). The parity transfor- 
mation and the associated minus sign are responsible for 
the identification {9 + tt, — n) = {9, n) already discussed 
in Section III A. 

To turn to a less trivial example, let us see how this 
works for the case of the tetrahedral phase. In this case 
r = T, the symmetry group of the tetrahedron. A tetra- 
hedron has 3 orthogonal 2- fold axes, and 4 3-fold axes. 
For the representative spinor given earlier 











(22) 



the z axis is aligned with one of the 2-fold axes. One 
easily verifies that z axis rotations through tt leave the 
spinor unchanged. Alternatively, we can align one of the 
3-fold axes with the z axis with the choice (see Fig. 3) 



/ \ 



(23) 



VV3/ 



(this is most easily seen by considering the Majorana 
polynomial, which has a root at z = 0). Now a rotation 
through ±2tt/3 is seen to reproduce the same spinor but 



T (the other non-trivial one-dimensional representation 
comes from changing the sense of the 3-fold axes) . The 
topological properties of the resulting space •^0(3)xC/(i) ^ 
and the implications for superfluid vortices in the tetra- 
hedral phase, were discussed in Ref. 17. 

For the rectangular phase F = D2 in general, but D4 
for the square case (77 = (2?! -I- 1) 7r/3). Here £)„ denotes 
the dihedral group. The only non-trivial phases occur in 
the latter case, as may be seen by considering the value 
r] = TT, when the four points lie on the equator of the 
Bloch sphere. Then we have 
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(24) 



and a ±7r/2 rotation about the z axis is seen to give rise 
to a —1. 

This construction generalizes readily to other ordered 
states of arbitrary spin once the corresponding stabilizer 
subgroups and phases are identified^^. Note that when F 
is discrete, as for the spin 2 phases other than the ferro- 
magnet, the order parameter manifold is 4 dimensional. 



2. Local geometry 

Next we turn to the local properties of the order pa- 
rameter manifold. The inner product naturally endows 
this space with a metric^^. Using the parametrization 
Eq. (20) consider two states (pufi and (pR'.e' related 



9' ^ 9 + d9 and R' ^ RR'f' with R^ 
an infinitesimal rotation corresponding to 



Sab — ^c^abc 



(^) 



We find the squared distance between these two states to 
be 

- = d9^ + ^aAgab - 2d9^p ■ (^S'^'Vo) 

(25) 

where the tensor gab is 



gab^l^'lisi'Ksi'^jcj,,. 



(26) 



with phase factors e 



±2m/3 



Eq. (25) makes it clear that if (f>QS^''''(f>Q ^ 0, rotations 
and phase changes are coupled together. There is an ar- 
bitrariness in the way the phase is apportioned between 
the two factors in Eq. (20) that amounts to a choice of 
gauge. This gauge structure was discussed in Ref. 11 for 
the case of the ferromagnet, the only one of the phases 
discussed in Section III with (^qS'^^^^q 7^ 0. The discus- 
sion of the general case is relegated to Appendix A. Our 



It is not hard to verify that main interest is in the other phases having (j>QS^'^\ 



0, 



these phases form a one-dimensional representation of for which Eq. (25) decouples into separate contributions 
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from the change in phase and the rotation, with the latter 
being characterized by the metric tensor gab- Note that 
gab is simply related to the nematicity Nab in Eq- (10). 
The notion of distance described by the metric tensor 
is left invariant i.e. preserved if states (jj^^ are mapped 



For the tetrahedral phase (x = ?7 = 7r/2) 3 = 2x1, 
showing that the order parameter manifold has a left and 
right invariant geometry. This is perhaps not surprising 
given the highly symmetric arrangement of points in the 
Majorana representation. For the rectangular phase (x — 
0), we see that in the case r] = 0, g — diag (3, 3, 0). The 
physical meaning is clear: because we have two points at 
either pole the stabilizer subgroup is 0(2), so that the 
order parameter manifold is only three dimensional. The 
same holds true for the spin 1 polar phase. 

We will see that the metric plays a crucial role in fixing 
the dynamics on the order parameter manifold. 



E. Conjugate variables 

Having characterized the order parameter manifold for 
the spinor condensates, we are almost ready to study the 
dynamics on that manifold. It remains to identify the 
conjugate variables. We expect these to be coupled by 
the first term of Eq. (4), which expresses the conjugacy 
of (j) and (j)^ . Substituting the parametrization Eq. (20) 
into that term gives 

= -dtd + 4iJfS^''>cl>^ 

(28) 

where ujt^a = \ (-^^'^t-^)cb components of the 

angular velocity. By analogy with rigid body dynam- 
ics, we refer to this as the 'body frame' angular veloc- 
ity: Eq. (28) shows that it is conjugate to the 'body 
frame' angular momentum (f)\S^^^ ip^, i.e. that of the im- 
rotated state. As is well known, the time integral of 
the term i(j)^dt4> alternative interpretation as the 

Berry phase associated with the time evolution of In 
this context the formula Eq. (28) appears in Ref. 28. 

Thus in the cases of interest where 4>\S^''^4'q = the 
variables conjugate to the rotations are nonzero only as 
one deviates from the order parameter manifold. To ac- 
count for these deviations we generalize the parametriza- 



by i? — > RR for some R, but not right invariant, under 
which g — s- RgR^ . 

Focusing now on the spin 2 case, we evaluate the metric 
for the state Eq. (17) 



(27) 
I 



N N 




S S 



FIG. 4: (color online) Distortion of the tetrahedral state 
shown in Fig. 2 due to boosts in the ^-direction 

tion Eq. (20) to 

cl>Bj,e^e''D^''\R)cl)i (29) 
where (f>i is a state with (j)\S^(l)i — I. Note that 

Rl is angular momentum in the 'lab frame'. To find the 
deviation corresponding to (pi , let us consider the state 

cj>i=U{l)B^'Kl)Po (30) 
where is some normalization factor, and 

Bi^\l)^c^^{^-{g-')jaSi'^y 

It is not hard to see that to quadratic order the normal- 
ization takes the form 

m) = ^-\{9-')ab^-lb+--- (31) 

and that to this order (j)\S^(j)i = I, as required. The effect 
of the Z-distortion on the tetrahedral state is shown in 
Fig. 4. 








2 + cos 77 — \/3 cos X sin V 

4sin2 a 
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Using the parametrization Eq. (29) in Eq. (28) gives 
then 

icj^'fdtcl) = -dtO + u> • /, 

to first order in I. 

We make the following aside. The matrix 
D^^^{R)B^^\l) that acts on (f>o is a polar decomposition 
of the Ds.o representation of an element of 5L(2,C)^^. 
Now SL(2, C) /Z2 is isomorphic to the connected Lorcntz 
group. This isomorphism has a beautiful physical coun- 
terpart. The Lorentz transformations have a natural ac- 
tion on the Celestial sphere, the space of light rays on the 
past (say) light cone. Further, the elements of SL{2, C) 



A- 



a b 
c d 



ad — be = 1, 



have a natural action on the Bloch sphere; 



A { ' ^ I , corresponding to a Mobius transformation on 



the stereographic coordinates z = a-f/aj. 

az + b 

Z 

cz + a 



(32) 



(the fact that A and —A correspond to the same Mobius 
transformation accounts for the Z2 quotient above). Re- 
markably, group elements related by the isomorphism 
mentioned above correspond to identical transformations 
of the sphere^°! The deformations of Eq. (30) arc just 
those generated on the Celestial sphere by a boost of 
the frame of reference, which in particular determines 
the aberration of the fixed stars. Thus the name 'stellar 
representation' sometimes used to describe the Majorana 
picture of spin states is more than picturesque. 

With the parametrization Eq. (29) we have accounted 
for six spin degrees of freedom (plus the superfluid 
phase): three rotation variables parametrizing the order 
parameter manifold and three conjugate variables. The 
remaining 4s — 6 variables required to specify the spin 
state do not correspond to any broken symmetries and 
must describe gapped modes. 

For the spin 2 case, the missing two degrees of freedom 
can be readily found by considering the ratio of the two 
(complex) 5Z/(2,C) invariants in Eqns. (18,19). 



/3 

J2 



(33) 



The powers are chosen so that the normalization in 
Eq. (30) drops out. By construction /z is unchanged for 
all I and R starting from the state Eq. (17), and could 
be expressed in terms of 77 and x- 

This concludes our discussion of the spin degrees of 
freedom in a spinor condensate. With this background, 
we will see that the derivation of the low energy La- 
grangian is extremely straightforward. 



IV. DYNAMICS NEAR THE ORDER 
PARAMETER MANIFOLD 

A. Low energy Lagrangian 

We are going to use the parametrization Eq. (29) in 

the Lagrangian Eq. (4). In doing so, we are treating the 
problem as as constrained dynamical system, in which de- 
viations from the ground state manifold associated with 
the gapped modes described in the previous section are 
assumed to be infinitely stiff. The Lagrangian for the 
spin degrees of freedom then takes the form 

Apin = U}fl - ^9abi^i,a^i,b " Hint(0 (34) 

where uj^^a = \^abc (-R.^9,i?)^^ i = x,y,z, in analogy to 
the earlier definition of uit^a- We could additionally allow 
for a variation 5p in the density of the gas, which is con- 
jugate to the phase 9, to describe the density modes, but 
will not do so here. In Eq. (34) we have not included the 
gradient of the conjugate variables (including 5p) in the 
part arising from the kinetic energy, as such terms can be 
neglected in the long wavelength limit of interest. The 
conjugate variables do however appear in the interaction 
term 



(35) 



The notation is of course chosen to emphasize the rigid 
body analogy. The precise form of the 'inertia tensor' 
will depend on the phase under consideration; we discuss 
the spin 2 phases for definiteness. In that case the in- 
teraction Hamiltonian has the form Eq. (7). Computing 
the quadratic variation of this expression with the conju- 
gate variables is facilitated by the SL{2, C) invariance of 
the third term: its variation is determined solely by the 
normalization in Eq. (31). We obtain 



ClSab 



tetrahedral phase 



ciSab ~ ^ {9 ^)ab rectangular phase 



(36) 



It is then straightforward to eliminate the I degrees of 
freedom using the equation of motion, obtained from 
Eq. (34) 



la — Iab'^t,b 



to obtain the final result 



■^spin — ^Iab<^t,al^t,b — 7:gab<^i,al^i,b (37) 



Eq. (37) represents the main conclusion of this work, be- 
ing the low energy Lagrangian for the spin degrees of 
freedom of the condensate. It takes the form of a sigma 
model in 3 -|- 1 dimensions with target space SO{3) (if we 
ignore the possibility of vortices in the superfluid phase 9 
the subtle global structure of the ground state manifold 
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discussed in Section III D 1 can be ignored) . In the case 
of the tetrahedral phase, the metric tensor g = 21L, and 
the Lagrangian Eq. (37) becomes that of the principal 
chiral model, having independent left and right SO{3) 
symmetries. 

An alternative form for Eq. (37) follows from noting 
that, if M = diag (mi, m2, ma) 



Mabi^H,a^ti,b = tr 



with 



M = (trM)l-2M 



Thus we have 
r ■ 

*^spin 



tr 



IdtR^dtR - gdiR^diR 



(38) 



By expressing the matrix elements of R in terms of an 
orthonormal triad Rab = (^6)a' ^ith • 66 = Sab this 
may be written 



1 ^ 

-Cspin = 2 X! {dteaf - 9a (VCa)' 



a=l 



(39) 



Recall that for the spin 1 polar phase and for the special 
value ?7 = in the spin 2 rectangular phase the metric 
tensor has one zero eigenvalue and two equal non-zero 
eigenvalues. As a result both g and / have fAuo vanishing 
eigenvalues, and Eq. (39) reduces to the usual 0(3)/0(2) 
sigma model. 



'^ti = df^tp. The linear equations of motion following 
from Eq. (41) are then 

dfltp - didigip = 0, 

a wave equation describing the propagation of three spin 
wave modes with velocities 



Va = 



a = 1,2,3 



For the spin 2 case Eq. (43) gives 



tetrahedral phase 
rectangular phase 



(43) 



with ga given by the diagonal elements of Eq. (27) with 
X = 0. For the square case (t/ = (2n + l)7r/3) we have 
g = diag (4, 1, 1) and 



vi = W4c: 



C2 

5 



^2 = ^3= WCi 



C2 
5 



(44) 
(45) 



These results check with Ref. 6, which also includes the 
normal phonon mode as well as the mode associated with 
variations of the r] and x parameters in Eq. (17). In the 
case of the rectangular phase this latter mode appears 
gapless in mean field theory, but as explained in Sec- 
tion III C 2 this is the result of an accidental degeneracy 
that does not persist in the next order of approximation. 



B. Equations of motion and spin wave spectrum 

We find the equations of motion corresponding to 
Eq. (37) by writing the variation 

{R^SR)^^^ = -IpcCabc 

This gives 

SiJf, = df^ip -ip xu}f, (40) 

Substitution into Eq. (37) leads to the equations of mo- 
tion. 

dt (/wt) - d, (gu^i) + wt X (luit) - i^z ^ (ff^,) = (41) 

It follows from their definition that u)^ satisfy the 
(Maurer-Cartan) equation 

d^ufi, - dt,u}^ + ^u>^ X = (42) 

Note that if a;^ is interpreted as a non-abelian gauge 
field, the above condition corresponds to vanishing field 
strength, and to the absence of topological defects. 

The equations of motion can be linearized by ignor- 
ing the right hand side of Eq. (42), allowing us to write 



V. DISCUSSION 

We have achieved our goal of providing a framework 
in which the parameters entering the low energy spin 
Lagrangian of an arbitrary ordered state of a spinor con- 
densate (with (p^S^'^^cpQ = 0) may be easily calculated. 
Though we focused on the spin 2 states, any other state 
can be treated by the same method once the problem of 
minimizing the mean field energy is solved. The exten- 
sion of the present formalism to spin ordered Mott insu- 
lating phases in which the phase variables are quantum 
disordered does not present any particular difficulties. 

Perhaps the most interesting problem that we have not 
addressed in detail relates to the character of topological 
defects in these systems. The occurrence of nonabelian 
stabilizer subgroups means that vortices have very novel 
characteristics^^'^^. We mention one consequence of our 
work for the quantum description of such vortices. The 
phase factors associated with elements of the stabilizer 
subgroups that were discussed in Section III D 1 will ap- 
pear in the path integral when vortices are present, as 
may be seen from Eq. (28). Consider an imaginary time 
path integral with fields obeying the boundary condition 

</.(r,r + /3) = <^(r,T) 
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If, as wc go from t t+i3, the field at a point r is subject 
to a rotation that evolves from R i?7, for 7 G F, the 
9 variable must increase 9-^0 — A(7) in order to ensure 
periodicity of the fields, leading to a phase factor e*''^^'''^ 
in the path integral {p is the density). Ref. 31 discusses 
the effec;t of these phases for the simplest case of the spin 
1 polar phase (or rather the Mott insulating phase based 
upon it), where they are ±1 and the defects are abelian. 
The nonabelian case remains unexplored. 
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APPENDIX A: METRIC ON THE ORDER 
PARAMETER MANIFOLD FOR (S) + 

In the general case the distance between states on the 
ground state manifold has the form Eq. (25). As men- 



tioned in the text, there is an arbitrariness in the appor- 
tioning of phase between the two factors of Eq. (20) that 
belies a natural gauge structure. In other words, Eq. (20) 
is unchanged under 

L'('')(i?) ^ D("'(i?)e^^^(^) 
9 9 + A{R) 

for some arbitrary function A{R). Let us define a (Berry) 
vector potential 

a = -i(j>lDUD(l}Q 

= • (SW)o (Al) 

where (• • • )o denotes an expectation in the state (j)^, 
'^a = \^abc {R^dR)^^, and d denotes the exterior deriva- 
tive (we find it convenient to use the language of differ- 
ential forms). The vector potential allows us to define 
a covariant derivative = d — ia. The metric tensor 
d(j}^ (g) dcf) can then be cast in the form 



iSid(f) = (f)ldD^ O dD^Q + d9 ® d9 + d9 ® a + a® d9 
= -(l)lD^DaD O D^V^D(f>o + {d9 + a) ^ {d9 + a) 

= uia<Si oJbQab + {d9 -h a) O {d9 + a) (A2) 

The two terms in Eq. (A2) are manifestly gauge invariant under the transformation Eq. (Al). The gauge invariant 
metric tensor Qab that appears in the first term takes the form^"" 

9ab = liiSi'^ - {Si%,si'^ - {Si%})o. (A3) 
In the second term, a by itself is gauge dependent, but the associated field strength is not 

da = -dw(S'("))o 

= ^ea6cWa Aa;6(S'('*))o (A4) 

where in the second step we have used the Maurer-Cartan equation Eq. (42). Eq. (A4) has a more familiar form, as 
may be seen by introducing a unit vector mo parallel to {S^^^)o. Then we have 
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= ea/37 {^^^^)aa' '^0,a' A (R^dR)^^, mo,/3'Wo,^ 

= Cafj^ {R^dR)^^, mo,a' A {R^dR)^^, mo,0' (R^R)^^, mo,Y 

— ^aff-fdRaa'nT'O.a' A dRpfj'TnQ^^' Rj-yiTTlQ^^' 

= eai3jdma A drriism-f (-^5) 
I 

where m = Rmo, and in the penultimate line we have ces is unity. As a result 
used the fact that the determinant of the rotation matri- ^ 

da = -eabc{S''a^)odmb A dvic, (A6) 
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which generahzes the Mermin-Ho relation in Eq. (2) to 
an arbitrary spin state. 

An important example is provided by the ferromagnet, 
for which the state (pQ is a fully spin polarized (coherent) 

state, and we have 

9 = ^ i^ab - mo^amofi) . 



Then the first term of Eq. (A2) takes the form 

u)a ® cObQab = ^dm (g) dm. 

The resulting metric sets the form of the Hamiltonian in 
Ref. 11. 
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